Viscous Withdrawal of Miscible Liquid Layers 
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In viscous withdrawal, a converging flow imposed in an upper layer of viscous liquid entrains 
liquid from a lower, stably stratified layer. Using the idea that a thin tendril is entrained by a 
local straining flow, we propose a scaling law for the volume flux of liquid entrained from miscible 
liquid layers. A long-wavelength model including only local information about the withdrawal flow 
is degenerate, with multiple tendril solutions for one withdrawal condition. Including information 
about the global geometry of the withdrawal flow removes the degeneracy while introducing only a 
logarithmic dependence on the global flow parameters into the scaling law. 
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Recent experiments on thermal convection with two 
layers of miscible liquids reveal several distinct regimes - 
an overturn regime with violent mixing of the layers, 
a doming regime where the interface undulates, and a 
stratified regime where the convection is largely stable 
with thin tendrils or sheets of one liquid entrained within 
the other Analogous steady-state entrained struc- 
tures arise in drainage flows 0, Q ) oil extraction Q , as 
well as viscous withdrawal of immiscible liquid layers, 
which occur in microfluidics [H, fiber coating 0] and en- 
capsulation of biological cells 7| . Recent works explor- 
ing the connections between thermodynamic phase tran- 
sitions and the topology transition that takes place at 
the onset of entrainment have noted that, in order for 
the entrained structure to be completely isolated from 
the large-scale flow dynamics, the shape of its base must 
be a power-law cusp [1, H, @, Intriguingly, experi- 
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12| on miscible entrainment also seem to show 
a robust cusp-like shape at the base of long-lived tendrils 
(see Fig. [1]). This suggests the entrained tendrils are iso- 
lated from the fluctuating, large-scale convection by the 
cusp-shaped base and are therefore able to remain sta- 
ble over many convection cycles. If true, this may even 
explain why hot-spots can persist over many convection 
cycles in the Earth's mantle [H, Motivated by these 
observations, we focus on the stratified regime in thermal 
convection of miscible layers and present a model that 
tests how the large-scale flow and topography anchor a 
thin cylindrical tendril. 

Because the large-scale flow is stabilized in the strat- 
ified regime, mixing between the layers is controlled by 
the volume flux of liquid entrained through the tendrils, 
Qq. Existing estimates of Qq assume that the veloc- 
ity field inside the tendril is uniforml y u pwards , flowing 
at a characteristic convection speed [111 [Hi [3l- How- 
ever, recent particle-image-velocimetry (PIV) measure- 
ments within the base of an anchored tendril reveal a 
stagnation-point velocity field, one more appropriately 
described by a characteristic strain rate E {s~^) instead 
of a characteristic velocity scale [lH . Viscous withdrawal 
experiments on immiscible layers suggest how an inte- 
rior stagnation point can arise 0. When the effect of 




FIG. 1: Axisymmctric tendrils entrained by a warmer, up- 
wards flow in a 2-layer thermal convection experiment with 
miscible oils. Convection is driven by heating the bottom of 
the tank and cooling the top. The lower fluid is less viscous 
than the upper and the spacing between tendrils is 0(10 cm). 
The dark vertical line on the left is a thermocouple [l^] . 

the entrainment penetrates deeply into the lower layer, a 
broad tendril forms with the interior moving uniformly 
upwards. This is the situation addressed by existing es- 
timates. When the effect of the entrainment penetrates 
weakly, a narrow tendril forms by drawing liquid inward 
within a thin layer below the interface, thus creating an 
interior stagnation point. Here we derive a new scaling 
law for the volume flux entrained, one which takes into 
account the interior stagnation point. 

The steady-state flow fields associated with non-linear 
interface deformations are difficult to obtain analytically. 
We therefore analyze the simplest case: the steady-state 
entrainment of a tendril from a deep lower layer by an 
axisymmetric viscous withdrawal flow imposed in an up- 
per layer (Fig. [5]). The entrained liquid is taken to be 
much less viscous than the exterior, so that there is only 
a weak feedback from entrained flow to the withdrawal 
flow. To ensure that the effect of entrainment pene- 
trates only weakly into the lower layer, we require that 
the layers are strongly stratified. With these simplifica- 
tions, the steady-state interface divides into three geo- 
metrically distinct regions: a long and slender entrained 
tendril where interior flow effects are dominant, a far- 
field interface where the hydrostatic pressure is domi- 
nant, and a transition region where the force balance 
changes smoothly from one form to another. 

First we will give a scaling argument that predicts the 
volume flux Qq. The system is best described by a cylin- 
drical coordinate system where z — corresponds to 
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FIG. 2: Sketch of surface R(z) and flow streamlines for an 
axisymmetric withdrawal flow in the upper liquid layer en- 
training a thin, cylindrical tendril from a deep lower layer. 

the height of the undisturbed interface and r = corre- 
sponds to the centerhne of the tendril. The steady-state 
tendril radius is defined as R{z). In the strongly strati- 
fied regime considered here, the stagnation point inside 
the tendril lies at a height Iz- The velocity field near 
the stagnation point is an axisymmetric straining flow, 
{—Er/2) Br + Ez Bz- To estimate Qo we note that, above 
H-zi the liquid from the lower layer is moving entirely up- 
wards so that Qo ~ {Ez)'kR^{z). Since it is possible to 
relate f^, the tendril radius near the stagnation point, 
to in a simple way described below, the volume flux is 

Oo - {Eiz)^e.\. 

Assuming the tendril is long and slender in the neigh- 
borhood of the stagnation point, the interior velocity field 
can be approximated via standard slender body type ar- 
guments [16]. The upwards contribution to the interior 
flow due to entrainment is then a uniform plug flow, Ez. 
The downward contribution is a pipe flow driven by hy- 
drostatic pressure Apg{R^{z) — r^)/(4^o), where Ap is 
the density difference between the two layers, g the grav- 
itational acceleration, and fiQ is the viscosity of the en- 
trained liquid. By definition, these two flows must cancel 
at the interior stagnation point {r = 0, z ^ £z) so that 



Eiz 



(1) 



It remains to determine iz- We identify £z as the am- 
plitude of the large-scale upwards deflection produced on 
the interface by the withdrawal flow. It is therefore set 
by a balance of /i-E, the upwards pull exerted by viscous 
stresses in the upper-layer withdrawal flow, where fi is 
the viscosity of the liquid in the upper layer, and Apg, 
the downwards pull of hydrostatic pressure. To simplify 
the notation for later analysis which uses £z as a charac- 
teristic lengthscale, we define 



4 = 2fiE/Apg. 
Relations ([T]) and ^ together yield 



\ Apg ) \ ApgJ- 



(2) 



(3) 



corresponding to the geometric mean of the viscous 
lengthscale associated with the upper layer flow and the, 
much smaller, viscous lengthscale associated with the less 
viscous, lower layer flow. The tendril "slenderness ratio" , 
^r/^z, assumed to be small at the beginning of our anal- 
ysis, is y/2fM)/p, which is indeed small for /io <C 

Finally, using ([T]) and we can rewrite Qo in terms 
of the local strain rate E and material parameters 



Qo = l67:coti^poE^/{Apgf 



(4) 



where cq is an 0(1) dimensionless entrainment coefficient 
which we determine later in the full analysis. This scaling 
law differs from existing estimates based on a character- 
istic convection velocity U 11 , [l^, [l3| ■ If we take S as 



the size of the convection cell, then ([4]) says Qo oc (U/S)'^ 
while existing estimates say Qo oc U^. We also empha- 
size that nothing in the above scaling argument depends 
on the withdrawal being axisymmetric. An analogous ar- 
gument for entrained 2-d sheets yields the same relations 
for £z (12) and the sheet thickness £r ^ , and the scaling 
law Qo/L ss E£z£r, for the volume flux per unit length. 

In deriving the scaling law for Qo O we essentially 
assumed the entrainment dynamics is controlled by a lo- 
cal straining flow. Since a straining flow has no inherent 
lengthscale, the dynamics is completely decoupled from 
the geometry of the large-scale withdrawal flow. In ther- 
mal convection, this feature would imply that the volume 
flux Qo has no dependence on the size of the large-scale 
convection cell. The last implication is counter-intuitive. 
After all, the tendril size at its base, where it joins on the 
interface, is controlled by how the interface levels out and 
becomes flat on the large lengthscales. This large-scale 
topography in turn depends on the specific geometry of 
the withdrawal flow, so physical intuition says that the 
entrainment dynamics, particularly Qo, should also de- 
pend on the global geometry. 

We next address this question by developing a long- 
wavelength model of the entrainment dynamics and find 
that a model which includes only the local straining flow 
is degenerate. For each value of Qo, there exists a con- 
tinuous family of tendril solutions, each with a different 
power-law shape at the base. The existence of a steady- 
state tendril thus does not require a specific power-law 
base shape. Including the large-scale withdrawal flow 
geometry removes this degeneracy. It also introduces a 
logarithmic dependence on the global dynamics into the 
scaling law for Qo (HI). 

Starting with the Navier-Stokes equations, and using 
standard slender-body approximations where effects pro- 
portional to £r/ £z are discarded at leading-order, we de- 
rive the following equation for the steady-state tendril 
shape 



Qo = {Ez)ttR^{z) 



nR^jz) dPp 
8p.o dz 



(5) 
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and interior pressure 



Pn 



2fiE (l 



dR 



R{z) dz 



Apgz. 



(6) 



In essence, equation ([5]) says that the unknown volume 
flux of liquid entrained into the tendril has two contribu- 
tions: a plug flow Ez induced by the withdrawal, and a 
pipe flow associated with non-uniform interior pressure. 
The expression for Pq ([6]) derives from the normal stress 
balance across the surface of the tendril. The second 
term is the familiar hydrostatic pressure difference. The 
first term is the large pressure necessary to keep a tendril 
from collapsing under the inward squeeze exerted by the 
exterior straining flow. Equations ([5]) and ([6]) together 
yield a second-order nonlinear ordinary differential equa- 
tion for the steady-state tendril shape R{z) and Qo- The 
derivation is analogous to the entrainment model for im- 
miscible liquids in [lo| . except that the surface tension 
contribution to Pq is replaced here by the hydrostatic 
pressure difference. 

The natural boundary conditions are conditions on 
how the tendril shape R{z) tapers to downstream as 
z —^ oo and how R{z) flares out into a flat interface up- 
stream as z — *■ 0. Both conditions correspond to asymp- 
totic balances of the governing equations. Far above the 
interface, the upwards plug flow Ez dominates and 



R{z) Roo{z) = ^Qo{co)/{ttEz) 



(7) 



where we have written Qo{co) to emphasize that in our 
analysis, we use the scaling law for Qo, so that cq 
is the only undetermined parameter. Near the tendril 
base, the effects of entrainment are negligible because 
the layers are strongly stratified. The tendril base shape 
is therefore determined by a balance of the hydrostatic 
pressure and the exterior viscous stress, thus 



R{z) -> Rs{z) 



z 



0. 



(8) 



Both the coefficient B and the exponent a are determined 
self-consistently with the full solution. 

Using boundary conditions (O, ([5]), and appropriate 
downstream structural stability modes, we numeri cally 
integrated equations ([5]) and ^ for tendril solutions [17 1 . 
We found that, for a fixed Qo, a continuous family of ten- 
dril solutions exist. Each solution's upstream shape is 
characterized by a different set of B, a values. Fig. [3] 
gives a few examples of the different tendrils possible 
for the same volume flux. The long-wavelength model 
is therefore clearly degenerate. 

Removing the degeneracy requires an additional con- 
dition on the upstream shape, one related to the large- 
scale topography. We next show this explicitly. In our 
toy model an axisymmetric withdrawal flow is generated 
in the upper layer by inserting a point force with strength 




FIG. 3: Rescaled tendril solutions at S/£z = 15, no/n = 0.1, 
Co = 3.1 with exponents a = 1.0, 0.5, 0.25, 0.052 (top to bot- 
tom). Inset shows tendril solution's power-law divergence. 

F at a height S above the undisturbed interface. This 
is clearly not a realistic withdrawal near the point force 
but, if S is large, it is a good approximation of real- 
istic withdrawal flows near the interface. The velocity 
fields in both fluid layers can be obtained by a variation 
of the method of images 18| and correspond to an ax- 
isymmetric straining flow with E = F / {2tt jiS"^) in the 
neighborhood of the base of the tendril. Since the effects 
of entrainment penetrate weakly into the lower layer, the 
large-scale interface shape should be well approximated 
by the deflection occurring when no liquid is entrained. 



R,{z) = S 7(34/2z)'/' 



1. 



(9) 



The fact that Ri{z) is controlled by S is the extra infor- 
mation needed to uniquely select a tendril solution. 

To incorporate this extra information, we require that 
at an unknown location Zs^ the tendril shape, as well as 
di?/dz and d^i?/dz^ equal the corresponding quantities 
in Ri{z) jl^. We solve for z^, a and B analytically as 
follows. Since Zg should be be near the base of the tendril, 
we approximate R{z) by the upstream shape Rs ([8]) and 
find that z^ « 0.609 4, a « 0.052, and B « 1.18 S. We 
can solve for the appropriate tendril solution numerically 
by varying the dimensionless entrainment coefficient cq so 
that R{z) and its derivative merge smoothly onto Ri{z). 
Fig. |4] shows an example of how cq is found numerically. 

It is interesting to know how the dimensionless entrain- 
ment coefficient cq depends on the global lengthscale, 
S/£z, and on the viscosity contrast /xo/m- We have var- 
ied S/i, from 10 to 10* and ^J.o/^J- from IQ-^ to lO"*^ 
and plotted the results in the inset of Fig. |4l The varia- 
tion of Co with these parameters is unexpectedly simple, 

well approximated by co = 2.3 log (^S/ {iz \/^o/Ai)) —0.74. 

The logarithmic dependencies of cq on S/£z and /io//^ are 
weaker than the power-law dependencies of Qo on £z and 
/io//i in the scaling law ([4]). This shows that including 
the global dynamics only weakly modifies the simple en- 
trainment law. 

The logarithmic dependence ultimately results from 
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FIG. 4: Rescaled tendril solutions R{z) for co = 3.4,3.1, 
and 2.9 (solid lines, top to bottom), together with the large- 
scale deflection solution Ki{z) (dotted) for Sjlz = 15 and 
Mo/m = 0.1. Only Co = 3.1 allows the two shapes to merge 
smoothly. Inset shows the values of co found from joining the 
numerical solutions smoothly onto Ri(z) in this way. 

the decaying exponential shape of the tendril at 
z ^ iz (HI). How this emerges is complicated because 
Co is essentially determined by demanding i?oo(-z) merge 
onto Rs{z) through a region where all the terms in the 
governing equation ([5|) are equally important. Our rough 
analysis shows that two requirements appear necessary 
for the observed Cq relation. First, the upstream ten- 
dril shape Rs{z) should assume a size consistent with 
Roo{z), the steady-state downstream shape. Second, the 
tapering of the tendril must be accomplished hy z ~ £z, 
with the sharp exponential decay in dRs{z)/dz appro- 
priately switching over to a gentler square-root decay in 
di?oo(2)/dz [13. 

In sum, we have now an improved scaling law for Qq, 
one which explicitly accounts for the dependencies on the 
global flow parameters 

where 71 =2.3 and 72 — —0.74. We have also analyzed 
tendril solutions associated with the withdrawal flows 
generated by a ring vortex and by a point sink and found 
the same qualitative outcome. The differences in geome- 
try between these withdrawal flows result in slightly dif- 
ferent tendril shapes and values of 71 and 72. Using (fTU)) 
we have estimated how quickly two miscible layers will 
mix in thermal convection and found results consistent 
with observed values. However, the strong dependence on 
E of pU)) makes a precise comparison difficult. To fully 
establish how the tendril persists over time, a complete 
numerical study as well as an experiment where the with- 
drawal is generated directly, e.g. by withdrawing liquid 
from a tube inserted into the upper layer 0, @, Q , instead 
of indirectly via thermal convection, are needed. 

Aside from providing an estimate for Qo that takes 
into account the interior stagnation point, our analysis 
also reveals that, although some information about the 
large-scale flow is necessary to specify a tendril shape, 



the dependence on the large-scale flow geometry is weak. 
This may be why tendrils in thermal convection experi- 
ments can remain stable over long periods of time despite 
fluctuations in the global convection which can alter the 
large-scale topography. 

The authors thank Anne Davaille, Sarah C. Case, Leo 
P. Kadanoff and Sidney R. Nagel for encouragement 
and helpful discussions. This work was supported by 
a GAANN Fellowship (L.E.S.) and NSF MRSEC DMR- 
0213745. 



[1] A. Davaille, J. Fluid Mech. 379, 223 (1999); M. Le Bars 
and A. Davaille, J. Fluid Mech. 499, 75 (2004). 

[2] S. Courrech du Pont and J. Eggers, Phys. Rev. Lett. 96, 
034501 (2006). 

[3] J. T. Jeong and H. K. Moffatt, J. Fluid Mech. 241, 1 
(1992). 

[4] S. Blake and G. N. Ivey, J. Volcanol. & Geotherm. Res. 
27, 153 (1986); J. R. Lister, J. Fluid Mech. 198, 231 
(1989). 

[5] A. M. Gafian Calvo, M. Perez-Saborid, J. M. Lopez- 
Herrera, and J. M. Gordillo, Eur. Phys. J. B 39, 131 
(2004); R. Suryo and O. A. Basaran, Phys. Fluids 18, 
082102 (2006); S. L. Anna, N. Bontourx, and H. A. Stone, 
Appl. Phys. Lett. 82, 364 (2003). 

[6] P. G. Simpkins and V. J. Kuck, Nature 403, 641 
(2000); J. Eggers, Phys. Rev. Lett. 86, 4290 (2001); 
E. Lorenceau, F. Restagno, and D. Quere, Phys. Rev. 
Lett. 90, 184501 (2003); E. Lorenceau, D. Quere, and 
J. Eggers, Phys. Rev. Lett. 93, 254501 (2004). 

[7] I. Cohen, H. Li, J. L. Hougland, M. Mrksich, and S. R. 
Nagel, Science 292, 265 (2001); J. L. Wyman, S. Kizilel, 
R. Skarbek, X. Zhao, M. Connors, W. S. DiUmore, W. L. 
Murphy, M. Mrksich, and S. R. Nagel, Small 3, 683 
(2007). 

[8] S. C. Case and S. R. Nagel, Phys. Rev. Lett. 98, 114501 
(2007). 

[9] I. Cohen and S. R. Nagel, Phys. Rev. Lett. 88, 074501 
(2002). 

[10] W. W. Zhang, Phys. Rev. Lett. 93, 184502 (2004). 
[11] A. M. Jellinek and M. Manga, Nature 418, 760 (2002). 
[12] A. Davaille, F. Girard, and M. Le Bars, Earth Planet. 

Sci. Lett. 203, 621 (2002). 
[13] A. M. Jellinek and M. Manga, Rev. Geophys. 42, RG3002 

(2004). 

[14] N. H. Sleep, Geophys. J. 95, 437 (1988). 

[15] A. Davaille, Poster, Focusing Stress in a Soft Interface 

Workshop, U. Chicago (2005). 
[16] A. Acrivos and T. S. Lo, J. Fluid Mech. 86, 641 (1978); 

G. I. Taylor, Proc. Int. Congr. Appl. Mech. Munich 11, 

191 (1964). 

[17] L. E. Schmidt and W. W. Zhang, unpublished. 
[18] S. H. Lee, R. S. Chadwick, and L. G. Leal, J. Fluid Mech. 
93, 705 (1979). 

[19] Here we have chosen the simpler option of matching Ri{z) 
and R{z) at a point instead of systematically matching 
them. Our numerics show the tendril merges smoothly 
onto Ri{z), so the errors appear small. The simpler pro- 
cedure also has the advantage of revealing parameter de- 
pendencies in an explicit way. 



